Spontaneous Emission Near Superconducting Bodies 
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In the present paper we study the spontaneous photon emission due to a magnetic spin-fiip 
transition of a two-level atom in the vicinity of a dielectric body like a normal conducting metal 
or a superconductor. For temperatures below the transition temperature Tc of a superconductor, 
the corresponding spin-flip lifetime is boosted by several orders of magnitude as compared to the 
case of a normal conducting body. Numerical results of an exact formulation are also compared to 
a previously derived approximative analytical expression for the spin-flip lifetime and we find an 
excellent agreement. We present results on how the spin-flip lifetime depends on the temperature T 
of a superconducting body as well as its thickness H . Finally, we study how non-magnetic impurities 
as well as possible Eliashberg strong-coupling effects influence the spin-flip rate. It is found that 
non-magnetic impurities as well as strong-coupling effects have no dramatic impact on the spin-flip 
lifetime. 
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It is well-known that the rate of spontaneous emission 
of atoms will be modified due to the presence of a dielec- 
tric body [Ij. In current investigations of atom micro- 
traps this issue is of fundamental importance since such 
decay processes have a direct bearing on the stability of 
e.g. atom chips. 

In magnetic microtrap experiments, cold atoms are 
trapped due to the presence of magnetic field gradients 
created e.g. by current carrying wires [^. Such mi- 
croscopic traps provide a powerful tool for the control 
and manipulatk)n of cold neutral atoms over microme- 
ter distances ||3|. Unfortunately, this proximity of the 
cold atoms to a dielectric body introduces additional de- 
cay channels. Most importantly, Johnson-noise currents 
in the material give rise to electromagnetic field fluctu- 
ations. For dielectric bodies at room temperature made 
of normal conducting metals, these fluctuations may be 
strong enough to deplete the quantum state of the atom 
and, hence, expel the atom from the magnetic microtrap 
Q • Reducing this disturbance from the surface is there- 
fore strongly desired. In order to achieve this, the use of 
superconducting dielectric bodies instead of normal con- 
ducting metals has been proposed Q . Some experimental 
work in this context has been done as well, e.g. by Nir- 
rengarten et al. 6] , where cold atoms were trapped near 
a superconducting surface. 

In the present article we will consider the spin-flip 
rate when the electrodynamic properties of the super- 
conducting body are described in terms of either a simple 
two-fluid model or in terms of the detailed microscopic 
Mattis-Bardeen [?} and Abrikosov-Gor'kov-Khalatnikov 
[§| theory of weak-coupUng BCS superconductors. In ad- 
dition, we will also study how non-magnetic impurities, 
as well as strong coupling effects accordiiig to the low- 
frequency limit of the Eliashberg theory [9|, will affect 
the spontaneous emission rate. 

Following Ref.[lo| we consider an atom in an initial 



state |i) and trapped at position = (0, 0, z) in vacuum 
near a dielectric body. The rate F^ of spontaneous and 
thermally stimulated magnetic spin-flip transition into a 
final state |/) is then 



Fs = ^0 



X Im[ V X V X G{rA,rA,uj) ]jk {n + 1) , (1) 

where we have introduced the dimensionless components 
Sj = {f\Sj/h\i) of the electron spin operators Sj with 
j = X, y, z. Here gs ~ 2 is the gyromagnetic factor of 
the electron, and G(r,r',Ci;) is the dyadic Green tensor 
of Maxwell's theory. Eq. ^ follows from a consistent 
quantum-mechanical treatment of electromagnetic radi- 
ation in the presence of an absorbing body [111 IT^ . In 
this theory a local response is assumed, i.e. the char- 
acteristic skin depth should be larger than the mean 
free path of the electric charge carriers of the absorbing 
body. Thermal excitations of the electromagnetic field 
modes are accounted for by the factor {n -\- 1), where 
n = i/(^g'i^/i^BT _ 2) g^j^fj ^ ^ 27ri^ is the angular fre- 
quency of the spin-flip transition. Here T is the tempera- 
ture of the dielectric body, which is assumed to be in ther- 
mal equilibrium with its surroundings. The dyadic Green 
tensor is the unique solution to the Helmholtz equation 



V X V X G(r, r', uj) - k^e{r, Lu)G{r, r', uj) = 6{r - r')l , 



(2) 



with appropriate boundary conditions. Here k = u/c is 
the wavenumber in vacuum, c is the speed of light and 
1 the unit dyad. The tensor G{r,r',uj) contains all rele- 
vant information about the geometry of the material and, 
through the relative electric permittivity e(r, w), about 
its dielectric properties. The fluctuation-dissipation the- 
orem is build into this theory 11 , 13] ■ 
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The decay rate F]^ of a magnetic spin-flip transition for 
an atom in free-space is well-known (see e.g. Refs.fl^). 
This free-space decay rate is Tg = TbS^, where = 
Mo (mb3s)^ / (Snh) and where we have introduced the 
dimensionless spin factor 5^ = + Sy + S^. The 
free-space lifetime corresponding to this magnetic spin- 
flip rate is = 1/F]^. In the present paper we only 
consider ^''Rb atoms that are initially pumped into the 
155*1/2, F — ^jTUp = 2) = |2,2) state, and assuming the 
rate-limiting transition |2, 2) —f |2, 1) in correspondence 
to recent experimen ts m , llil - [13, • The spin factor is 
52 = 1/8 (c.f. Ref.[3) and the frequency is = 560 
kHz. The numerical value of the free-space lifetime then 
is = 1.14 X 10^5 s. 

In the following we will consider a geometry where an 
atom is trapped at a distance z away from a dielectric slab 
with thickness H. Vacuum is on both sides of the slab, 
i.e. e{r,uj) = 1 for any position r outside the body. The 
slab can be e.g. a superconductor or a normal conducting 
metal, described by a dielectric function e(w). The total 
transition rate for magnetic spontaneous emission 



(3) 



can then be decomposed into a free part and a part purely 
due to the presence of the slab. The latter contribution 
for an arbitrary spin orientation is then given by 
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with the atom-spin orientation dependent integrals 
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Cm (x) 
Cm (ix) S> , (6) 



where the scattering coefficients are given by [If 

2 _ gix H/z 

Cn{x) = rp{x) 
Cm{x) = rs{x) 



1 - r2(a;)e*^^/^ ' 

2 _ ^ixH/z 

1 - r2(a;)e^^^/^ ' 



(7) 



(8) 



The electromagnetic field polarization dependent Fresnel 
coefficients are 



rp{x) 



e{uj)x - ^ {2kzf{e{w) ~ I) + 



e{uj) X + ^ {2kz)'^{ €{uj) - 1) + 



, (9) 



X - V(2fcz)2(eH-l) + x^ 

rs(x) = , . (10) 

X + ^{2kz)^{e{u;)~l) + x^ 



For the special case H = 00, the integrals in Egs.® 
and ([6]) are simply a convenient re- writing of Eqs.(8)- 
(12) in Ref.[17]. Note that « 2/|| provided kz < 1. 
Throughout this article, we use the same spin-orientation 
as in Refs.|^[ll], i.e. 



5*2 = 52 and Sx = 0. 
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FIG. 1: tb of a trapped atom near a superconducting film as 
a function of the temperature T/Tc- The solid as well as the 
dashed-dotted line correspond to the two-fluid mode and the 
Gorter-Casimir temperature dependence. We use Al(0) = 
35 nm and S{Tc) ~ 150 ^m [23|, corresponding to niobium. 
The critical temperature is Tc = 8.31 K For T/T^ > 1 

we put a2{T) ~ but ai{T) = 2/uj^j.o5(T^f ■ The dashed line 
corresponds to a film made of gold described by the dielectric 
function given by Eq. (|15p . The upper most graph (dotted 
line) shows to the lifetime T%/{n + 1), i.e. the free-space 
lifetime with r| = 1.114 x 10^^ 

As the total current density responds linearly and lo- 
cally to the electric field, the dielectric function can be 
written 



e{Lj) = 1 h I 



(11) 



Here a{T) = cri(r) + iu2{T) is the complex optical 
conductivity. We may now parameterize this complex 
conductivity in terms of the London penetration length 
Al(T) = ,/l/w^oo-2(r) and the skin depth 5{T) = 
\/ 2 / Ld noai (T). In this case, the dielectric function is 
e(w) = 1 - l/fc2Ai(T) + i2/PS^{T). If, in addition, 
we consider a non-zero and sufficiently small frequency 
in the range < ut <^ ojg = 2A(0)/?i, where A(0) is 
the energy gap of the superconductor at zero tempera- 
ture, the current density may be described in terms of 
a two- fluid model [31 . The London penetration length 
is Ai(T) = yJns{T)/nQ and the skin depth is 

^(r) = 5{Tc)l yJnn(T)/no. Here the electron density 
in the superconducting and normal state are ns{T) and 
n„(T), respectively, such that ns(T) + n„(T) = tiq and 
ns(0) = n„(T > Tc) = no [l9| . A convenient summary 
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of the two-fluid model is expressed by the relations 



luniT) 

no 



a^m^aJ^, (12) 



no 



where ct„ = ai{Tc) and ctl = l/a;/ioA|^(0). Considering, 
in particular, the Gorter-Casimir temperature depen- 
dence [2^ for the current densities, the electron density in 
the normal state is n„(T)/no = (T/Tc)^. For niobium we 
use S{Tc) = y^2/w^oO'ti ~ 150 /xm as cr„ « 2x lO'^(Slm)^^ 
and Ai,(0) — 35 nm according to Refs.(2l|. In passing, we 
remark that the value of (t„ as obtained in Ref. is two 
orders of magnitude larger than the corresponding value 
inferred from the data presented in Refs.[21||. 

The lifetime tb = l/r^ for spontaneous emission as 
a function of T is shown in Fig. [1] for 7? = 0.9 /xm (solid 
line). We confirm the observation in Ref.fl^ that for 
temperatures below Tc and for H = 00 (dash-dotted 
line), the spin-flip lifetime is boosted by several orders of 
magnitude. In Ref. [ist . the spin-flip lifetime was, how- 
ever, calculated by making use of the approximative and 
analytical expression 



^i(r) 



1 



3/2 



(T) (kzY 



, (13) 



valid provided Al(T) < d{T) and Al(T) < z < A. Com- 
paring this analytical expression with the numerical re- 
sults as presented in Fig. [U based on the exact equations 
Eqs.([51)- pII)) . we find an excellent agreement. This obser- 
vation remains true when ai (T) and (T2 (T) are obtained 
from more detailed and microscopic considerations to be 
discussed below. For temperatures T/Tc > 1 we can 
neglect the (T2(T) dependence and, for S{T) ^ z, we 
confirm the result of Ref.fS], i.e. 



2eow 



(14) 



For T ~ Tc we have to resort to numerical investigations. 

In contrast to the traditional Drude model, more real- 
istic descriptions of a normal conducting metal in terms 
of a permittivity include a significant real contribution to 
the dielectric function in addition to an imaginary part. 
One such description is discussed in Ref.[2§], where 



e{uj,T) = 1 



-I- I 



(15) 



and nv{T) = Q.OMl{T/ef J^^'^ dx x^e'/ie' - if eV us- 
ing a Bloch-Griineisen approximation. Here 9 = 175 K 
for gold. The plasma frequency is fiLOp — 9 eV. For tem- 
peratures T ~ 0.25Tc, we observe that Eq. p5|) leads to 
(Ti(r) ~ (72 (T), and that for lower temperatures (72(T) 
will be the dominant contribution to the conductivity. 
For temperatures T/Tc > 1, in the use of Eg. p^ we can 



set (72 (T) ~ when calculating the lifetime. For a bulk 
material of gold this leads to almost two orders of mag- 
nitude longer lifetime as compared to niobium since for 
gold S{Tc) ~ 1 /im, using the parameters corresponding 
to FiglU This finding is in accordance with Eg. p^ . As 
seen from FiglU for a thin film and for T/Tc > 1 we 
find the opposite and remarkable result, i.e. a decrease 
in conductivity can lead to a larger lifetime. 

A much more detailed and often used description of 
the electrodynamic properties of superconductors than 
the simple two-fluid model was developed by Mattis- 
Bardeen 0], and independently by Abrikosov-Gor'kov- 
Khalatnikov [1], based on the weak-couphng BCS theory 
of superconductors. In the clean limit, i.e. I 3> 5oi where 
I is the electron mean free path and ^0 is the coherence 
length of a pure material, the complex conductivity, nor- 
malized to cr„ = (Ji(Tc), can be expressed in the form 

n 

cr(r) f°° dx , fx + nuj\ , , 

= / T- tanh g{x) 



where g{x) 



A(T) 
(x2 - 



dx 

- — tanh 

flLO 



A^(r) 



2fcsT 



(16) 



hu}x)/uiU2 and ui = 
Here, the 



A2(T), U2 = .yix + huj)^ - A2(T) 
well-known BCS temperature dependence for the super- 
conducting energy gap A(T) is given by [2^ 



In 



v/(?i^i3)2 + A2(0) /A(0) 



dx 



^x^ + A^T) 



tanh 



^x^ + A^T) 
2kBT 



(17) 



where lud is the Debye frequency and A(0) — 
3.53 kBTc/2. For niobium, the Dcbyc frcquency_is Hlod = 
25 meV. According to a theorem of Anderson |23, [1^ , the 
presence of non-magnetic impurities, which we only con- 
sider in the present paper, will not modify the supercon- 
ducting energy gap as given by Eq. p7|) . The complex 
conductivity will, however, in general be modifled due to 
the presence of such impurities. 

In the dirty limit where I ^0, the complex conduc- 
tivity has been examined within the framework of the mi- 
croscopic BCS theory (see e.g. Ref.fs^). In this case, the 
complex conductivity, now normalized to ctl, can conve- 
niently be written in the form 



a(T) 



dx 

— tanh 
2 



flUJ 



1 



2kBT 



1 



U2 — ui + iTi/t U2 -I- Ml — ih/r 



A(T) 



dx 

— tanh 
2 

g{x) + 1 



2kBT 



9{x) - 1 



U2 — ui + ih/r U2 + ui + ifi/r ^ 



(18) 
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Here we choose r such that h/r A{0) = tt^o/^ = 13.61, 
corresponding to the experimental coherence length = 
39 nm and the mean free path 1{T ~ 9K) = 9 nm. The 
normaUzation constant is ctl = 1.85 x 10^"^ (f^m)"^ cor- 
responding to Al(0) — 35 /im for niobium |2l| . 



values of tb as obtained using the two-fluid model for 
-ff = oo as presented in Fig. [1] and the corresponding 
result as shown in Fig. [3] shows, for our set of physical 
parameters, that the two-fluid model overestimates tb 
with three order of magnitude. 
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FIG. 2: The complex conductivity a{T) = ai{T) + ia2{T) as a 
function of the temperature T/Tc with R/tA(0) = 13.61 
The solid line corresponds to Eq. (I18|) . and the dashed-dotted 
line corresponds to Eq. (US} and a2{T) = 0.25 x 10"V2(r). 

As the temperature decreases below Tc, Cooper pairs 
will be created. Despite a very small fraction of Cooper 
pair for temperatures just below T^, the imaginary part of 
the conductivity as given by Eq. (|18p exhibits a vast in- 
crease (c.f. Fig. [2]). Furthermore, due to the modification 
of the quasi-particle dispersion in the superconducting 
state, there is an increase in ai{T) as well just below Tc- 
This is the well-known coherence Hebel-Schlichter peak 
[2^ . In contrast to the simple Gorter-Casimir tempera- 
ture dependence, both Egs. fTB]) and (|18p describe well the 
presence of the Hebel-Schlichter peak, with a peak height 
less then 8ct„ for both cases (c.f. Fig.[2|), at least for the 
values of the physical parameters under consideration in 
the present paper. In the opposite temperature limit, i.e. 
T < Tc, numerical studies of Eq. HH) show that <Ji{T) 
decreases exponentially fast. As seen in Fig. [2l the imag- 
inary part of the conductivity, on the other hand, is more 
or less constant for such temperatures. 

In passing we observe that there is only a minor dif- 
ference in (J2{T) as obtained from Eqs. (fT6|l and (|18p re- 
spectively. For temperatures around the peak value of 
the Hebel-Schlichter peak, cri(T) obtained from Eq. (fT8|) 
is, however, approximatively twenty percent larger than 
(Ti (T) as obtained from Eq.([T6]). This difference has, nev- 
ertheless a small effect on the lifetime tb- Hence, com- 
puting Tb using Eqs.pB]) or (fTS]) for the complex con- 
ductivity, we realize that the presence of non-magnetic 
impurities have no dramatic impact on the lifetime for 
spontaneous emission (see Fig. A comparison of the 
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FIG. 3: Tb of a trapped atom near a superconducting bulk as 
a function of the temperature T /Tc. The other relevant pa- 
rameters are the same as in Fig. [T] The solid line shows the 
lifetime tb using the microscopic BCS theory, i.e. Eq. (|18|l . 
The bold dotted line corresponds to the Mattis-Bardeen the- 
ory, i.e. using Eq. (|16|) . The dashed-dotted line shows the life- 
time Tb using Eq. (O, with ft/rA(0) = 13.61. The dashed 
line corresponds to a film made of gold described by the di- 
electric function given by Eq. (I15|) . The upper most graph 
(dotted line) shows to the lifetime T%/{n + l), i.e. the free- 
space lifetime with r£ = 1.114 x lO^''. 

For finite values of the lifetime r and for non-magnetic 
impurities we can also investigate the validity of the two- 
fluid model approximation in terms of the lifetime tb 
for spontaneous emission processes. As we now will see, 
there are large deviations between the microscopic theory 
and the two-fluid model approximation, in particular for 
small temperatures. According to Abrikosov and Gor'kov 
(for an excellent account see e.g. Ref.[29] and references 
cited therein), the density of superconducting electrons 
is given by 

provided that rA(0)/?i <C 1. We can now compute 
the dielectric function Eq. pT|) using Eq. p^ . We find 
that (72 (T) I (Jl obtained in this way agrees well the cor- 
responding quantity obtained from Eq. (|18p . There is, 
however, a considerable discrepancy between the two- 
fluid expression for cri (T) / ul and the corresponding ex- 
pressions obtained from the microscopic theory as given 
by Eq. (fT8|) . The numerical results for the lifetime in this 
case are illustrated in the dashed-dotted lines in Fig. [3] 
and in Fig. U) 
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FIG. 4: Tb of a trapped atom near a superconducting film as 
a function of the temperature T/Tc with H = 0.9 /^m. The 
other relevant parameters and labels are the same as in Fig. [S] 



Since we are considering low frequencies < co <^ ujg = 
2A(0) /h, strong coupling effects can now be estimated by 
making use of tlie low-frequency limit of the Eliashberg 
theory Q and its relation to the BCS theory (see e.g. 
Ref. 3l|). The so-called mass-renormalization factor Zjv, 
which in general is both frequency and temperature de- 
pendent, is then replaced by its zero-temperature limit, 
which for niobium has the value ~ 2.1 |31||. Using the 
strong-coupling expressions for the optical conductivity 
in a suitable form as e.g. given in Ref. we then 

find that the complex conductivity cr(r)/cr„ is rescaled 
by cr„ (7n/ZN with the lifetime of non-magnetic impu- 
rities rescaled by r — > t/Z]\[. The change in the lifetime 
for spontaneous emission can then e.g. be inferred from 
the relation Ea. (|13p . and we find only a minor decrease 
of Tb by the numerical factor « 0.69, which also 

agrees well with more precise numerical evaluations. 

The lifetime for spontaneous emission exhibits a mini- 
mum with respect to variation of the thickness H of the 
superconducting film. This fact is illustrated in Fig. [S] 
Below the minimum at Hmin ~ 0.1 /.tm, a decrease of the 
thickness H leads to an increase of lifetime in proportion 
to . This happens despite the growth in polarization 
noise because the region generating the noise is becoming 
thinner as it is limited by iJ, and not Xl{T). Eventu- 
ally, the lifetime reaches the free-space lifetime is as 
H tends to zero. On the other hand, for large H, i.e. 
H ^ ^{T), the lifetime is constant with respect to H, 
giving the same result as for an infinite thick slab. In the 
region between, i.e. Xl{T) ^ H < S{T), the lifetime is 
proportional to H. Numerical studies show that a non- 
zero a2 (T) is important for a well pronounced minimum 
of Tb as a function of H. 

Some experimental work has been done using a su- 
perconducting body, e.g. Nirrengarten et al. Here 
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FIG. 5: Tb as a function of the thickness H of the film. 
The complex conductivity is computed applying Eq. (|18p with 
ft/TA(0) = 13.61. Other relevant parameters are the same 
as in the solid line in Fig. [1] In the limit H = 0, i.e. no 
slab at all, the lifetime is simply Tg/in + 1) = 7.34 • 10^^ 
s for the parameters under consideration. The skin depth 
5{Tc/2) — 80.5 fj,m is indicated by a dashed-dotted line, and 
the London length Al(Tc/2) = 89.2 nm is indicated by a 
dashed line. 



cold atoms were trapped near a superconducting sur- 
face. At the distance of 440 /zm from the chip sur- 
face, the trap lifetime reaches 115 s at low atomic den- 
sities and with a temperature 40 /xK of the chip. We 
believe the vast discrepancy between this experimental 
value and our theoretical calculations must rely on ef- 
fects that we have not taken into account in our analysis. 
The use of a thin superconducting film may lead to the 
presence vortex motion and pinning effects in (see e.g. 
Refs. 32, d^l). The presence of vortices will in general 
modify the dielectric properties of the dielectric body. 
If we, as an example, consider a vortex system in the 
liquid phase in a finite slab geometry, one expects a 
strongly temperature dependent ai (T) with a peak value 



ai{T) ~ 1.3 X lOy H'^[fim]iy[kH z]nm '32\. Close to this 
peak o'i(T) ~ a2{T), and for v ~ 560 kHz we find a life- 
time for spontaneous emission two orders of magnitude 
larger than a film made out of gold with the same ge- 
ometry. It is an interesting possibility that spontaneous 
emission processes close to thin superconducting films 
could be used for an experimental study of the physics 
of vortex condensation. This possibility has also been 
noticed in a related consideration, which has appeared 
during the preparation of the present work [3^ . There 
are also fabrication issues concerning the Nb-0 chem- 
istry (ssj which may have an influence on the lifetime for 
spontaneous emission. 

To summarize, we have studied the rate for sponta- 
neous photon emission, due to a magnetic spin-flip tran- 
sition, of a two-level atom in the vicinity of a normal 
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conducting metal or a superconductor. Our results con- 
firms the conclusion in Ref.[l^, namely that the corre- 
sponding magnetic spin-flip lifetime will be boosted by 
several orders of magnitude by replacing a normal con- 
ducting film with a superconducting body. This con- 
clusion holds when describing the electromagnetic prop- 
erties of the superconductivity in terms of a simple two- 
fluid model as well as in terms of a more detailed and pre- 
cise microscopic Mattis-Bardeen and Abrikosov-Gor'kov- 
Khalatnikov theory. For the set of physical parameters 
as used in Ref.fll] it so happens, more or less by chance, 
that the two-fluid model results agree well with the re- 
sults from the microscopic BCS theory. We have, how- 
ever, seen that even though the two-fluid model gives 
a qualitatively correct physical picture for spontaneous 
photon emission, it, nevertheless, leads to large quantita- 
tive deviations when compared to a detailed microscopic 
treatment. We therefore have to resort to the microscopic 
Mattis-Bardeen Q and Abrikosov-Gor'kov-Khalatnikov 
theory in order to obtain precise predictions. We 
have also show that non-magnetic impurities as well as 
strong-coupling effects have no dramatic impact on the 
rate for spontaneous photon emission. Vortex conden- 
sation in thin superconducting films may, however, be of 
great importance. Finally, we stress the close relation be- 
tween the spin-flip rate for spontaneous emission and the 
complex conductivity, which indicates a new method to 
experimentally study the electrodynamical properties of 
a superconductor or a normal conducting metal. In such 
a context the parameter dependence for a bulk material 
as given by Eq. lfTB]) may be useful. 
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